Energy funneling in a bent chain of Morse oscillators with long— range coupling 
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A bent chain of coupled Morse oscillators with long-range dispersive interaction is considered. 
Moving localized excitations may be trapped in the bending region. Thus chain geometry acts like 
an impurity. An energy funneling effect is observed in the case of random initial conditions. 
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1. INTRODUCTION 

Nonlinear excitations (solitons, discrete breathers, in- 
trinsic localized modes, etc.) have been drawing increas- 
ing attention over recent years and are widely believed 
to be responsible for several effects in molecular chains, 
such as charge and thermal conductivity, energy transfer 
and localization, etc. (see reviews in 0,0, IS EM El e.g.). 

Initially, the geometrical features of the polymers and 
biopolymers were essentially neglected and energy trans- 
fer and localization was mostly attributed to inhomo- 
geneities and impurities H, [ j , l a , ITfl HH IT^ | or nonlin- 
ear excitations 0, El IT3. Tl5[1lrl Il7j| . Also, discreteness 
plays an important role for the localization of these exci- 
tations. The inhomogeneities have been modelled by dif- 
ferent masses at various chain sites 

U El EOLky chan ges 
in the coupling between molecular sites 0, 13 or by dif- 
ferent on-site potentials E3 as we U as conformational 
defects |l7j . In general, impurities have been shown to 
act as filters governing the progression of incoming exci- 
tations. Thus both reflection, trapping and transmission 
of incident moving discrete breathers through the impu- 
rity region can occur H, El El 0, El • Similar effects 
have been observed through collisions between moving 
discrete breathers, thus [8llT^. ITr| showed how station- 
ary large amplitude discrete breathers, on the average, 
absorb energy from colliding breathers of smaller ampli- 
tude. Thus the large amplitude breather may play a role 
similar to that of an impurity |12| . 

Recently, both long-range dipole-dipole interaction 
helicity and curvature gHHIll have 

been included in the nonlinear transport theory, as well 
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as combinations of these effects [H US IHIH El It has 
been shown that chain geometry induces effects similar 
to those of impurities pa I2I I2I |25 | . 

Special attention has been paid to models of biologi- 
cal macromolecules , such as proteins |S El and DNA 
IS 113 IH El EE H3 These are obvious choices for 
more complex geometric models, as their structure is cru- 
cial for their functionality |3lL l32| . A widely used model 
was presented by Peyrard and Bishop [33j in the con- 
text of statistical mechanics. In biological environments, 
thermal fluctuations are always present and have been 
considered in |!J Hfl yj, |3f| |3g , e.g. In these refs. it was 
shown that solitons or discrete breathers can be gener- 
ated from initial random thermal fluctuations. 

The aim of the present work is to study the interplay 
between chain geometry and long-range interaction in 
an augmented Peyrard-Bishop model. We show how a 
new mechanism for energy accumulation in the system 
— funneling — may be provided by the geometry of the 
chain. We consider a simple approximative description 
of the long-range intersite coupling by modelling it as a 
dipole-dipole-like interaction. Using an attractive long- 
range interaction, we study the effect of the geometry of 
the bent chain on the dynamics. The particular shape 
of the bend turns out to make no qualitative difference 
in terms of trapping and funneling. We therefore choose 
a simple wedge-shaped geometry. As initial conditions 
we use discrete breathers as well as randomly distributed 
fluctuations [1113. 

In Sec. we introduce the model, which includes 
the long-range interactions and the chain geometry. In 
Sec. |21 we investigate breather dynamics in the system 
and in Sec. random initial conditions are considered. 
Sec. El summarizes our results and contains a discussion. 
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2. THE MODEL 

We consider a one dimensional lattice of Morse oscil- 
lators with the Hamiltonian density 

rLn = + -n \ U n ~ u n-l) 

+ [exp(-w n ) - 1] - - J nm u n u m , (1) 

m 

where prime indicates m ^ n in the summation. The 
Hamiltonian for the system becomes H = Yln^-N^m 
where the total number of sites is Nt = 2N +1. In 
Eq. (JTJ the first term is the kinetic energy at the n'th site. 
Then follows a harmonic potential interaction between 
neighboring sites, C being the dispersion parameter. An 
on-site Morse potential (shown in Fig. ^| describes the 
atomic interaction. Finally, there is a summation of long- 
range interactions in which the coefficients are given by 



where denotes the fixed wedge angle (see Fig. 0). We 
note that the geometry of the chain only comes into play 
through the long-range interactions. In fact, earlier stud- 
ies of the long-range effect in curved molecular chains 
show that the exact form of the additional dispersion is 




FIG. 2: Wedge chain with opening angle, 6. 
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(2) 



where Jo is a strength parameter, and r„ denotes the po- 
sition of the n'th site. In our model the distance between 
neighboring sites, |r n+ i — r„|, is constant and normal- 
ized to unity. 
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FIG. 1: The Morse potential, V(« n ) = [exp(-lt„) - l] 2 . 

From the Hamiltonian we get the equations of motion 
tin + C (2u n - u n _i - U n+ i) 



-2e-«™ [e""" - 1] - JnmUm = 



(3) 



At the ends of the molecule we use the free boundary 
conditions 

M-JV-l = U-jv, ( 4 ) 
un+i = u N . (5) 

The wedge-shaped chain is given by 



r n = {x n ,Vn) = (nsin-,|n|cos-), 



Throughout the pap er we use C — 0.075 (a typical 
value for DNA [iM 12313) an d Jo — 0.5 (based on esti- 
mates of the dipole moment H^). 

Fig. |3 gives a detailed picture of the long-range inter- 
action by plotting the interaction coefficient J nm versus 
n, for fixed (a) m = —3 and (b) m — —2. We see that 
the closer m gets to the bend at n = 0, the higher the 
shoulder in the J„ m -profile becomes. This feature sug- 
gests an analogy in which the bend acts as an impurity. 
Away from the bend this effect rapidly drops off. 
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FIG. 3: Long-range interaction coefficients Jnm for fixed (a) 
m = -3, (b) m = -2. 6 = 35°. 

The long-range interaction in curved chains may be 
represented as follows: 
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where the first summations on the right hand side corre- 
spond to the inhomogeneous dispersion (seen in Fig EJ) , 
while the second summation, in which V® s = — J2' m Jnm 
is introduced, corresponds to an effective on-site poten- 
tial [23. The potential V® s has the double-well profile, 
shown in Fig. ^ 
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FIG. 4: Effective potential, V n Eff 
number, n. 
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In the following we shall see how the impurity, or effec- 
tive on-site potential, introduced by the bend can reflect, 
trap or transmit incoming excitations. 



3. ENERGY TRAPPING 

In this Section we consider the interaction between the 
bend and the incoming localized excitations. 

A 4th order Runge-Kutta solver is used to simulate 
Eq. (JJJ on a chain with Nt — 301 sites. A stepsize in 
time of 0.005 ensures conservation of the Hamiltonian 
to a relative accuracy of 10 -10 throughout. We use a 
Gaussian initial condition 



Un(t) = Aexp — k((n — v) — vtY 



(6) 



where site v denotes the initial position of the center of 
mass. 

In the following we use the velocity v = 0.2, the width 
k = 0.2 and the amplitude A = 0.5, because they turn 
out to provide the right balance between nonlinearity and 
dispersion to allow the initial condition to evolve rapidly 
into a discrete moving breather. Insertion of Eq. © 
with these parameter values into the total Hamiltonian 
gives H — 0.13. After some initial radiation, the moving 
breather turns out to possess the energy, H s w 0.08. 

In the following sections we present the numerical sim- 
ulations of the chain dynamics. 



3.1. Breather dynamics 



only a slight decrease of the velocity after passage of the 
center region is observed. In contrast for a stronger bend, 
9 — 95°, a considerable part of the excitation is trapped 
at the tip of the wedge n = 0. 




FIG. 5: Contour plots for the evolution of the Hamiltonian 
density, Tin, for iV T = 301, J = 0.5, v = 0.2, k = 0.2, 
A = 0.5 and v = —100, yielding H — 0.13. 5 equidistant 
contours H = 0.005 . . . 0.05. (a) 6 = 140° (transmission), (b) 
9 = 95° (trapping). 

Very strong bends (smaller wedge angles) turn out to 
result in reflection of the incident breathers. Such scat- 
tering properties and their dependence on the strength of 
the bend-induced impurity, are similar to those of a lin- 
ear impurity |37| and those of large amplitude breathers 
acting as an effective impurity [lj|- We stress that the 
specific shape of the bend does not affect the scatter- 
ing properties of the bent chain. Thus similar properties 
were observed in a parabolic chain |38| . 

To analyze the processes in detail, we calculate the 
central energy, H c , in 21 sites around n = (21 being 
a typical span of the denaturation bubble of the DNA 
molecule [EHH) 
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(7) 



where 7i n is given by Eq. (JTJ- The results are shown 
in Fig. [S] In the transmission case with a small bend, 
Fig. EJ a )> nearly all the energy leaves the region. In 
the trapping case with a stronger bend, Fig. EJb), the 
trapped energy is stabilized at about H c /H s w 64%. 
Thus only part of the energy is trapped. 
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FIG. 6: Relative central energy, H c /H a , versus time for the 
simulations in Fig. El (a) 9 = 140°, (b) 9 = 95°. 



In Fig. |S1 we show contour plots for the evolution of 
the Hamiltonian density. A weak bend with 9 — 140° 
has no noticeable effect on the breather (Fig. [SJi) and 



The trapped energy portion, H c /H s , calculated at time 
= 800 is shown in Fig. [7| as a function of the wedge an- 
gle, 9. For smaller and larger 0-values, the energy is lost 
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through reflection and transmission, respectively. Effi- 
cient trapping is found for intermediate wedge angles, 
90° < 9 < 107°. The optimal wedge angle for trapping 
is seen to be around 9 = 95°. 




Wedge angle, 



FIG. 7: Relative central energy, H c /H s , at t — 800 versus 
wedge angle, 0. System parameters as in Fig.|S] 



4. MULTIPLE BREATHER DYNAMICS 

In Fig. [S] we show the trapping of multiple breathers. 
A first Gaussian pulse (7) is launched at site v = —100 
at t = 0. At t = 800 both the displacements and the 
velocities are set to zero, w„(800) = and it„(800) = 0, 
outside the bent region, \n\ > 15, to remove radiation. 
This "cleaned" chain is now used as an initial condition 
for a new simulation, in which we add a second identi- 
cal Gaussian pulse (77) launched at site v — —87. Like 
in other systems the interaction between two breathers, 
or a breather and an impurity, depends strongly on the 
relative phase. We choose v = —87 for the launching of 
this second pulse to obtain maximal trapping. Using the 
same procedure, a third identical Gaussian pulse (777) 
is launched at t = 1600, now at v = —89. As seen in 
Fig. |S] we essentially succeed in trapping 3 breathers at 
the tip of the wedge chain. Some energy transmission 
is observed when breathers 7 and 77 are trapped, while 
reflection occurs at the trapping of breather 777. By 
successively handling the initial conditions as described 
above, we avoid radiation which, when reflected at the 
boundaries, distorts the numerical simulations. 

The corresponding energy evolution for the central 
sites is shown in Fig. EI The ability of the system to trap 
energy at the bending region is evident, even though more 
radiation is observed as the number of trapped breathers 
increases. As noted also in connection with Fig. EJb), 
the first incident breather, 7, loses about 36% of the to- 
tal energy before trapping. For the following breathers, 
77 and 777, both of the corresponding losses amount to 
50%. Thus the possibility for trapping more energy at the 
chain bend by additional incoming breathers may seem 
exhausted due to an effective saturation. 
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FIG. 8: Trapping of breathers 7, 77 and 777 at n = 0. 
Contour plot for the evolution of the Hamiltonian density, 
Tin, with 5 equidistant lines H = 0.005 . . . 0.05. N T = 301, 
Jo = 0.5, v = 0.2, k = 0.2, A = 0.5. 7 : v = -100 at t = 0, 
77 : v = -87 at t = 800, III : v = -89 t = 1600. 
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FIG. 9: Central energy, 77 c , versus time corresponding to 
Fig. El 



5. FUNNELING 

The trapping of breathers observed in the previous Sec- 
tion suggests that the bend may funnel energy from the 
surrounding region. In order to study this in detail, we 
now explore the dynamics of the chain in the case of 
random initial conditions. 500 different realizations with 
zero displacement, u n — 0, and velocities normally dis- 
tributed with zero mean, <?i„>= and standard devia- 
tion <ju n = 0.17 are used, corresponding to a Hamiltonian 
77 « 1.41. 

Random initial disturbances may create nonlinear 
localized excitations, interacting with each other and 
with the effective inhomogeneity caused by the bend 
Hi OH [3 EM Hll- In 0, the number of generated solitons 
were found to depend on the temperature of the system, 
T, by the power 1/3. Here we find that collision of the 
nonlinear excitations may result in exponential growth 
of the oscillation amplitude at the collision site. We ob- 
serve this phenomenon in Fig. 1101 in a bend chain. Here, 
a sudden increase of the amplitude of the center site, uq, 
occurs after about 130 time units of bounded oscillations. 
This unbounded growth of amplitude implies energy lo- 
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calization at the center site. 
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FIG. 10: Displacement at center site, Mo, versus time on bent 
chain with Nt = 99 and 9 = 95°. Initial displacements: 
u n = for all n. Random initial velocities generated with 
normal distribution: <u n >= and au = 0.17. H=1.41. 
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Site, n 

FIG. 12: Wedge chain, 8 = 95° , Nt = 99. Occurrence of 
amplitudes above threshold, u n > 10, versus site n, until 
max. t — 10.000. Initial displacements: u n = for all n. 500 
random initial velocities with normal distribution: <ii n >— 
and au n = 0.17, H=1.41. 



Our results for a straight chain are depicted in Fig. ITU 
showing a histogram of the occurrences of sites with a 
displacement above the threshold value u n = 10, cor- 
responding to the value for DNA opening used in |39| . 
The simulations were discontinued when this threshold 
value was exceeded. If this event did not occur within 
10.000 time units, a "no-occurrence" was registered. 
"No-occurrence" happened in 24 out of 500 simulations. 
The threshold transgressions is seen to be uniformly dis- 
tributed along the chain. 

On a wedge chain with bending angle 6 = 95°, identi- 
cal initial conditions gives the remarkably different result 
shown in Fig. ^| Here, 93 (out of 491) threshold trans- 
gressions occurs in the center region —1 < n < 1. Only 
9 "no-occurrences" were registered. Thus energy local- 
ization, implied by unbounded growth of amplitude, is 
observed in the vicinity of the tip of the wedge which 
therefore acts as an energy funnel. 
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Site, n 

FIG. 11: Straight chain, 9 = 180°, N T = 99. Occurrence 
of amplitudes above threshold, u n > 10, versus site n, until 
max. t = 10.000. Initial displacements: u n = for all n. 500 
random initial velocities with normal distribution: <u„>— 
and au n = 0.17, H=1.41. 

For lower standard deviations, cr^ n , correspond- 
ing more realistically to physiological temperatures, a 



smaller number of simulations produce threshold trans- 
gressions, within 10.000 time units. Thus, for a^ n — 
0.116, H « 0.7, corresponding to the temperature T k 
37° C, only about 1% of the 500 simulations led to thresh- 
old transgression. The few transgressions that did occur, 
were found in the vicinity of the tip of the wedge. 



6. CONCLUSION 

On a bent chain of Morse oscillators we find that mov- 
ing discrete breathers may be trapped at a bending point 
in the presence of dipole-dipole-like longe-range interac- 
tion. Thus the role of the geometry for the dynamics is 
analogous to that of an inhomogeneity. At the bending 
point, several incident discrete breathers may be trapped. 
However, energy is lost to radiation and a saturation ef- 
fect seems to limit the total trapped energy in the vicinity 
of a given bending point. 

For random initial conditions modelling thermal fluc- 
tuations, the tendency to unbounded amplitude growth 
in the vicinity of the bending point is substantially am- 
plified. Thus energy localization is implied in this region 
which therefore acts as an energy funnel. 

The use of a nonlinear potential is crucial for obtaining 
the energy funnncling effect in our model, since no ampli- 
tude growth is observed in a linear approximation. The 
plateau of the characteristic Morse potential allows for 
the breaking of the hydrogen bonds, e.g., in the molecule 
to be modelled. In contrast, a linear approximation with 
a parabolic potential produces too powerful an attraction 
for this effect to take place and would therefore not be 
an adequate description of chemical bonds. 

We also note that the attracting nature of the long- 
range interaction is crucial for the occurrence of ampli- 
tude growth. Ongoing work on this effect uses a more 
accurate dipole-dipole interaction term and includes the 
effect of the twisting of the dipoles occurring along the 
strands of DNA molecules, e.g. 
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